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ABSTRACT 

The radio luminosity function (RLF) of radio galaxies and radio-loud quasars is often mod- 
elled as a broken power-law. The break luminosity is close to the dividing line between the 
two Fanaroff-Riley (FR) morphological classes for the large-scale radio structure of these ob- 
jects. We use an analytical model for the luminosity and size evolution of FRII-type objects 
together with a simple prescription for FRI-type sources to construct the RLF. We postulate 
that all sources start out with an FRII-type morphology. Weaker jets subsequently disrupt 
within the quasi-constant density cores of their host galaxies and develop turbulent lobes of 
type FRI. With this model we recover the slopes of the power laws and the break luminosity 
of the RLF determined from observations. The rate at which AGN with jets of jet power Q 
appear in the universe is found to be proportional to Q~ 16 . The model also roughly predicts 
the distribution of the radio lobe sizes for FRII-type objects, if the radio luminosity of the tur- 
bulent jets drops significantly at the point of disruption. We show that our model is consistent 
with recent ideas of two distinct accretion modes in jet-producing AGN, if radiative efficiency 
of the accretion process is correlated with jet power. 

Key words: galaxies: jets - galaxies: active - galaxies: luminosity function - radio contin- 
uum: galaxies 



1 INTRODUCTION 

The large scale structure, or lobes, produced by the jets of radio 
galaxies and radio-loud quasars shows a large range of different 
morphologies (for an overview see Muxlow & Garrington, 19911. 
Nevertheless they can be grouped into two main classes according 
to their surface brightness structure ( Fanaroff & Riley] [1974) . Ob- 
jects in Fanaroff-Riley class I, FRI for short, are brightest at their 
centres while FRII-type objects are edge-brightened. The two FR 
classes are also broadly separated in terms of their radio luminos- 
ity. The original dividing line was placed at 5 x 10 25 WHz -1 at an 
observing frequency of 178 MHz, but was later shown to depend 
on the properties of the host galaxies (Ledlow & Owen] [1996). 

Higher angular resolution observations revealed that the more 
luminous FRII-type sources contain laminar jets extending from 
the AGN to very bright radio hotspots enclosed in diffuse, low sur- 
face brightness lobes. The morphology is consistent with the inter- 
pretation of the hotspots as the impact sites of the jets on the ambi- 
ent gas and the lobes as the 'waste baskets' for the material and en- 
ergy transported by the jets (Scheuer 1974). The FRI class is made 
up of sources with different morphologies. However, the majority 
shows jets that appear laminar in the innermost region inflating tur- 
bulent lobes after passing through a flare point. A minority of FRI 
sources shows a 'fat double' morphology which is reminiscent of 
the diffuse lobes of the FRII class, but without hotspots ( |Owen &] 
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|Laing|[7989[|Owen & White] [199T) . In this paper we concentrate 
on the turbulent FRI-types. 

The luminosity dividing the FR classes is remarkably close to 
the break in the radio luminosity function (RLF) (e.g. |Dunlop &] 
|Peacock[|1990) . While this has been noted many times, most work 
on the RLF has focussed on its cosmological evolution and its sep- 
aration, not into FR classes, but into contributions from starburst 
galaxies, flat and steep spectrum objects (lackson & Wall. [1999) . 
Both FR classes are part of the steep spectrum population which 
also forms the parent population for the relativistically beamed flat 
spectrum objects (flat-spectrum quasars and BL Lacs). Advances in 
the modelling of the dynamics (e.g. Falle 1991 , Kaiser & Alexan- 
der^! 997] hereafter KA) and the synchrotron emission from the ra- 
dio lobes of individual FRII-type sources (Kaiser et al. 1997] here- 
after KDA; IBlundell et al.ll 19991 IManolakou & Kirk |2002) now 
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allow us to investigate the shape of the RLF above the break lu- 
minosity in detail. In this paper we present such an analysis. We 
argue that FRI-type sources evolve out of FRII-type objects with 
weak jets. This allows us to construct a self-consistent explanation 
for how the shape of the RLF arises from the combination of both 
FR classes and the luminosity evolution of individual sources. 

Differences in the properties of radio sources were also shown 



to be present beyond the radio waveband (e.g. |Hine & L ongair 
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ticular the difference in the emission properties of the AGN itself 
have let to the idea of different accretion modes in low and high 
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luminosity AGN {Falcke et ak"l |2004; Hardcastle et al. 2007}. It is 
important to note that there is no one-to-one correspondence be- 
tween the emission properties and postulated accretion mode of the 
AGN itself and the FR radio morphology. However, we show that 
our model developed here extends the idea of different accretion 
modes to include the FR dichotomy, if AGN with higher nuclear 
luminosities also give rise to more powerful jets. 

The observed distribution of the sizes of radio lobes of FRII- 
type objects is difficult to reconcile with evolutionary models. 
There appear to be too many small sources with sizes of tens of pc 
to a few kpc that cannot all evolve into objects with lobes tens and 
hundreds of kpc across (O'Dea & Baum 1997 1. We show that our 
model can explain these observations, since many small sources 
starting out with an FRII-type morphology 'drop out' of the ob- 
served samples by developing turbulent jets before their lobes grow 
to large sizes. 

In section [2] we review the luminosity and size evolution of 
FRII-type objects. We present a simple model for the transition to 
and subsequent evolution of FRI-type objects in section|5] Section]?] 
contains the construction of the RLF using these combined models. 
The distribution of lobe sizes for the FRII-type objects is discussed 
in section [5] In section [6] we fit our model into the observational 
constraints from wavebands other than the radio and discuss the 
proposal of two distinct accretion modes for jet-producing AGN. 
Section[7]provides a summary of our main findings. Throughout we 
use the definition of the spectral index a of synchrotron emission 
as F v oc v~ a . 



2 EVOLUTIONARY PHASES OF A FRII-TYPE 
RADIO-LOUD AGN 

In this section we investigate the evolution of the radio luminosity 
of radio-loud AGN of type FRII as a function of source age and 
lobe size. This evolution is governed by the energy loss processes 
of the radiating relativistic electrons in the radio lobes. We are us- 
ing the models of KA and KDA for this analysis and the relevant 
mathematical expressions are derived in the appendix. 

The key relation is given in equation (A19) . Unfortunately, 
a general analytic solution of the integral in this equation is not 
possible. However, we can simplify the analysis considerably by 
setting the exponent of the initial power-law energy distribution of 
the relativistic electrons, m, to 2. This choice implies that the spec- 
tral index of the radio emission is fixed to oc = (m— l)/2= 0.5 
in the case of negligible radiative energy losses of the relativistic 
electrons. We will show below that the spectral index also cannot 
become steeper than a — 1 while the lobes remain observable in 
this idealised model. While these restrictions clearly limit the ap- 
plicability of the model, the reduction in mathematical complexity 
allows us some important insights in the luminosity evolution. 

Setting m = 2 we solve the integral in equation l |A19[ l, 



„1+E 
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We further simplify the analysis by assuming that the equation of 
state of the lobe material is dominated by the relativistic electrons 
and the magnetic field, i.e. k = and T/ = 4/3. Hence the exponent 
e = a\/3 (see equation |A2TJ. 

The value of x m j n determines the fraction of the lobe volume 
that contributes to the emission at the observing frequency v. As we 
show in the appendix | A3 1 x m ^ n is governed by the energy losses of 
the relativistic electrons in the lobe. The energy losses are due to the 



adiabatic expansion of the lobe, the emission of radio synchrotron 
radiation and the inverse Compton scattering of CMB photons. It is 
not possible to find an explicit expression for x m i„ taking into ac- 
count all loss processes. However, while adiabatic losses affect the 
electrons at all times, the radiative loss processes usually dominate 
at different times during the source evolution. Hence we can define 
different phases during the evolution of a radio-loud AGN of type 
FRII during which different analytical expressions for their radio 
luminosity can be derived. 



2.1 Adiabatic losses 

We shall see below that energy losses due to the emission of syn- 
chrotron radiation affect the source evolution mainly at early times. 
Inverse Compton losses dominate at late times. Therefore in be- 
tween these two regimes a source can go through a phase where 
only adiabatic losses are important. This is the simplest possibility 
and we therefore consider it first. 

Neglecting radiative energy losses amounts to setting g(x,t) = 
0. Equation l |A26fr can then be solved to give 

Ymax / 

Unless the observing frequency is very high, we expect y v <C y ff 
and so the luminosity of the source in this phase is 



(2) 



L v = 
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With equations \KL\ and \PsA\ we can express the source age t as 
a function of the lobe length and the pressure inside the lobe. The 
luminosity is proportional to p 7 / 4 D 3 . For our self-similar model the 
volume of the lobe, V, is proportional to Z) 3 and the energy density 
of the magnetic field inside the lobe is proportional to p. Therefore 



7/4 

we recover the well-known result that Ly °= uJ V (e.g. Longair 



1994) for the minimum energy condition of a uniform magnetisec 
plasma inside V. 

Note also that the luminosity evolution in the adiabatic regime 
can be positive for p < 8/7. Sources located in atmospheres with 
a comparatively flat density distribution increase in luminosity as 
they grow, if radiative energy losses can be neglected. The spectral 
index in this regime is a = 0.5. 



2.2 Inverse Compton losses 

The mathematical description of the inverse Compton losses of the 
relativistic electrons is the same as for synchrotron radiation losses 
with the energy density of the magnetic field, uq, replaced with the 
energy density of the CMB photon field, mcmb ( see equation |A22[ l. 
While mb decreases with time, mcmb remains virtually constant 
over the lifetime of the jet flow. It is therefore inevitable that the 
inverse Compton losses will eventually dominate over synchrotron 
losses. In this regime we can neglect uq compared to mcmb and 



therefore equation (TA231 reduces to 
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This equation still takes into account adiabatic energy losses. Sub- 
stituting this into equation {A26} and solving for x m ; n yields 

1/04 
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:t 3(P-8)/[4(5-P)] > 



(6) 



The time-dependent part of this expression is given by l/y v ? 
From equations l |A4) and < |A14[ > we find that 

1 

For any reasonable choice of p, the function 1 /y v t decreases with 
increasing /. Hence at late times in the source evolution we can use 
a binomial expansion to simplify equation {5} and find the luminos- 
ity in this regime as 

3m e c 2 /« 



14\/3mcmbV 



(7) 



The luminosity in this regime always decreases as the source grows. 
Also note that the spectrum of the source now has a spectral index 
of a = 1 . 



Model parameter 


Value 


r, 


4/3 


r\ 


5/3 


P 


0,2 


A 


4 


Ymin 


1 


fmax 


10 6 


m 


2 


Q 


10 38 W 


P 


lCT^kgnr 3 


a 


2 kpc 


T 


10 7 K 


"CMB 


4x l(T 14 Jm- : 



Table 1. Parameters of the fiducial model 



2.3 Synchrotron losses 

As pointed out in the previous section, the energy density of the 
magnetic field, uq, decreases as the source grows. So at early times 
in the source evolution will dominate over mcmb and the lu- 
minosity evolution is governed by the synchrotron losses. In the 
expression for g (x, t) we now neglect mcmb an d so we arrive at 

l/o 3 



= 1- 



3a2tn e c 1 
4(J T JvtUB 
for a\ 7^ 3/5 and 
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for a\ = 3/5. For adiabatic and synchrotron losses the time- 
dependent part of the expression forx m ; n is given by 1/Yv/mb- From 
equations \A4\ , \Al l\ and l |A14[ l we now find 

-J—oefCMVWS-P)]. (10) 
Yv/t'B 

For P > 8/7 the importance of synchrotron losses decreases as the 
source grows, while for p < 8/7 the opposite applies. In environ- 
ments with a steep density distribution the pressure, and therefore 
the strength of the magnetic field, inside the lobes decreases fast 
enough so that the continued injection of freshly accelerated rel- 
ativistic electrons overcompensates the synchrotron losses. In less 
stratified atmospheres the synchrotron losses are so severe that the 
addition of new relativistic electrons cannot compete with them. 

Analogous to the discussion in the section above, we can de- 
rive an expression for the radio luminosity when the second term in 
the bracket in equation <j8j is small. This limit corresponds to young 
sources for (3 > 8/7 and to older sources for p < 8/7. In both cases 
we can use another binomial expansion and get 

m e c 2 f„ 
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Q = constant. 



(11) 



We arrive at the same result for the case a\ =3/5 by expanding the 
exponential function to first order. For p > 8/7 sources start out at a 
constant luminosity which is independent of the density parameter 
pa$, i.e. it is independent of the properties of the environment it 
is located in. As the source grows, synchrotron losses become less 
important, the luminosity starts to decrease and the source enters 
the regime where only adiabatic losses are important. For p < 8/7 
sources start in the adiabatic regime and increase in luminosity until 
synchrotron losses dominate and the luminosity takes on a constant 
value, which is again independent of the source environment. 



Whenever the luminosity is dominated by synchrotron losses, 
the spectral index is again given by a = 1 . 



2.4 Overall source evolution 

In general the gas density in the environments of radio-loud AGN 
is not well fitted with a single power-law. A P-model of the form 



(12) 



[l + (r/a) : 
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usually fits the X-ray emission from the hot gas in elliptical galax- 
ies, galaxy groups and galaxy clusters (e.g. Fukazawa et al.[|2004| >. 
The lobes of a radio-loud AGN with FRII-type morphology will 
therefore first expand in an atmosphere with roughly constant den- 
sity until it reaches a size comparable to a. At this point the source 
evolution changes, since for larger sizes the density of the external 
medium is approximated by a power-law as given in equation \Al) . 
The evolution of a source in a changing density profile is described 
in detail in Alexander ( 2000). Here we are interested mainly in the 
overall luminosity evolution of sources as they grow as we will use 
this information in constructing the radio luminosity function and 
the distribution of lobe lengths. 

It is useful to calculate order-of-magnitude numerical esti- 
mates for some of the quantities predicted by the model. For this 
purpose we set model parameters to physically meaningful, fidu- 
cial values. The model input parameters are summarised in table^ 
while model quantities derived from these inputs are collected in 
tableU 

In our model all sources start out with an FRII-type morphol- 
ogy. They first propagate through the central regions of their en- 
vironments where the density is essentially constant. During this 
phase 1, p = and the lobe length increases as (see equation|A2|> 



D 



(13) 



The luminosity evolution of the source in phase 1 is initially 
governed by adiabatic losses only which implies an increasing lu- 
minosity proportional to £) 2 / 3 (see equationBJI. Synchrotron losses 
become more important though as the source grows and the lumi- 
nosity levels out. For the values of our fiducial model we find from 
equation l|9j that x m j n ~ 0.5 for D = a. Synchrotron losses there- 
fore halve the luminosity of the source compared to the case of 
pure adiabatic energy losses. 

For very small sources we would need to take into account 
synchrotron self-absorption as demonstrated by the spectral shape 
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Density regime 


Parameter 


Value 


independent 


r 


3/4 




Pa 


i .4 x iu J m 




fn 


1.2 x 10 I2 s 2 kg- 1 m- 2 


(5 = 2 


CI 


1.5 




"\ 


3/2 




a 3 


-3/2 




®4 


1/2 




fp 


0.11 




Jy 


'■ z lGHz/ j m 




A 


3.4xlO- 17 (^)" I/2 J 1 ^m 1 /4 s 2 kg -i 


p = o 


C\ 


1.7 




"I 


3/5 




a-i 







as, 


4/5 




fp 


0.11 




fy 


7.3( 5fe ) 1 / 2 j'/4 m -3/4 




A 


3.4xl0-' 7 ( n | E )- 1/2 J 1 / 4 m 1 /Vkg- 1 



Table 2. Model parameters derived for the fiducial model in the density 
regimes discussed in the text. Note that the value of the derived parameters 
often does not depend strongly on the value of p. 



of Giga-Hertz Peaked Spectrum (GPS) sources. However, we are 
mostly concerned with larger sources, at least comparable to the 
sizes of Compact Steep Spectrum (CSS) objects, and so we neglect 
self-absorption. 

Note that the radio spectrum in phase 1 is predicted to steepen 
as the source grows and the luminosity tends to a constant value. 
Our model therefore predicts that young sources should have rela- 
tively steep spectra as confirmed for CSS objects. The steep spec- 
tra arise because of the growing importance of synchrotron losses 
in the comparatively dense core region of the source environments. 
However, the CSS objects are not 'frustrated' by an exceptionally 
high density. The uniform density of the environment simply limits 
the expansion speed in such a way as to allow synchrotron losses 
to accumulate. Even sources in less dense environments will go 
through this phase. It is the shape of the density distribution rather 
than the absolute value of the density itself that leads to a steep 
radio spectrum for CSS objects. 

During phase 1 sources with weaker jets may become turbu- 
lent and develop an FRI-type morphology. This is discussed below. 
Sources retaining laminar jets and their large-scale morphology of 
type FRII grow eventually larger than the core radius a and enter 
phase 2. The external density distribution is now approximated by 
a power law. The lobes grow faster than before, but the exact re- 
lation between the lobe length D and the source age t depends on 
the value of the power law exponent (3 (see equation [A2]l. In phase 
2 synchrotron losses become less important as the lobes grow pro- 
vided that (3 > 8/7. The spectrum flattens again and the luminosity 
evolution follows equation {3}. 

The last phase, phase 3, is reached when the energy density of 
the magnetic field in the lobe falls below the energy density of the 
CMB. The lobe continues to expand in the same way as in phase 2, 
but the luminosity evolution steepens according to equation l|7j. At 
the same time the spectrum also steepens again. 

The end of the source lifetime is reached when the jet flow 
stops altogether. Depending on the fuelling mechanism for the 
AGN activity and the jet production mechanism this may occur 
during any of the evolutionary phases discussed above. However, 
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Figure 1. Sketch of the luminosity (solid line) and spectral index (dashed 
line) evolution of a lobe with an FRII-type morphology. The transitions 
from one regime to the other are of course expected to be less sharp than 
shown here. 



for most powerful jets sustaining an FRII-type morphology a typ- 
ical jet lifetime of 10 8 years allows the source to evolve through 
all three phases. A schematic representation of the luminosity and 
spectral index evolution of the lobe as a function of its length is 
presented in Figure[T| 



3 FRI SOURCES 

For the purposes of this paper, we only consider FRI-type sources 
with turbulent jets. We do not consider FRI-type sources with a 'fat 
double' morphology (e.g. Owen & Laing 1989; Owen & White 

[T99"IT >. 

3.1 Transition to an FRI-type morphology for sources with 
weak jets 

The development of an FRII-type structure, i.e. a laminar jet flow 
inside a lobe, depends crucially on the pressure inside the lobe ex- 
ceeding the pressure in the source environment. The jet propagates 
through the low density lobe and its turbulent disruption by fluid 
instabilities is unlikely. When the pressure inside the lobe becomes 
comparable to the ambient pressure, the lobe surface becomes un- 
stable and the denser gas initially replaced by the lobe expansion 
starts to refill the lobe volume (e.g. Alexander, 2002J. For a roughly 
isothermal external atmosphere with temperature T the pressure 
distribution has the same profile as the density distribution, i.e. 

PO 



Px 



[l + (r/a) 2 ] 



lP/2' 



(14) 



with po = kftTp/m, where m is the mean mass of the gas particles. 
The pressure is higher closer to the centre of the distribution and 
so the replacement of the lobe material starts here after pressure 
equilibrium is reached. This gives rise to the central 'pinching' of 
the lobe observed in some FRII-type sources (Baldwin 1982). 

For sources with lobes of length D < a, the lobe pressure 
decreases rapidly, p oc Z)~ 4 / 3 , while the pressure of the ambient 
medium is roughly constant. Once the lobe extends beyond a, the 
pressure of the source environment decreases as while the lobe 
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pressure is proportional to £)(~ 4 ~P)/ 3 . Since for the great majority 
of sources 1 < (3 < 2, the ratio of ambient and lobe pressure remains 
almost constant. In other words, once a lobe escapes the core region 
of its environment without reaching pressure equilibrium, it will re- 
main overpressured for the rest of its evolution. 

The dense external gas replacing the lobe will eventually reach 
the jet flow itself. At this point the protective function of the lobe 
disappears and the jet can develop instabilities on its surface lead- 
ing to its turbulent disruption. At this point the source morphology 
changes from FRII to FRI. The limiting case of a lobe just escap- 
ing from the core region and therefore avoiding the disruption of its 
embedded jet is given by a lobe with a radius D/A = a at the point 
in time where p = py. Using equation \A4\ we find the limiting jet 
power 

Q min =(^) 3/2 A^ T VW- (15) 
\ m fpJ 

For our fiducial model we derive 2 m ; n = 3.3 x 10 37 W, close to 
the dividing line between the two FR classes found by Rawlings & 
|Saunders]jl991) from the observed correlation of radio luminosity 
and the strength of optical line emission from the AGN itself. Jets 
with energy transport rates below <2 m ; n will start out with an FRII- 
type morphology, but turn into FRI-type objects once their jets are 
disrupted. More powerful jets in the same environment retain their 
FRII-type morphology throughout their lifetime. 

Both the jet power and the properties of the ambient den- 
sity distribution determine whether a source makes the transition 
to an FRI-type morphology or not (see also Bicknell 1994 , Gopal- 
|Krishna & Wutal |2TM) |Alexander| [20001 on this point). We there- 
fore expect that sources with jet powers close to Q m \ n and located 
in environments with an asymmetric density distribution may de- 
velop a hybrid morphology. The fact that only few such objects are 
known, all of which have radio luminosities close to the demarca- 
tion line between the FR classes ( Gopal-Krishna & Wiita , [2000| >, 
argues for mainly symmetrical source environments. 

3.2 Dynamics and luminosity of an FRI source 

Turbulent FRI-type jets often appear to emerge as laminar flows 
from their AGN. They then develop a flare point where the jet sud- 
denly widens. The flare point is usually the brightest feature of the 
large-scale radio structure in these objects. After passing through 
the flare point, the jet flow appears turbulent and resembles smoke 
rising out of a chimney (e.g. Muxlow & Garrington, 1991 1. It seems 
natural that the jet decelerates strongly in the flare point and that the 
synchrotron emitting particles are mainly accelerated here (Baum| 
|et~al| [T9971 |Laing & Bridie] |2002b|a[ [2504) l. Jet deceleration and 
particle acceleration in the flare point strongly suggest an interpre- 
tation of this point as a strong shock within the jet flow. 

Models that relate the radio luminosity of the jet to the tur- 
bulent jet flow after passing through the flare point are difficult to 
construct and mainly rely on conservation equations rather than de- 
tailed dynamics (Bicknell 1994, 1995). It is also possible to infer 
details of the jet flow from high resolution radio observations of 
turbulent jets ( |Laing & Bridle||2002bja||2004| >, but a fully analytic 
model for the luminosity of FRI-type sources based on the jet dy- 
namics equivalent to the model for FRII sources described above 
and in the appendix is not available. Despite this lack of a detailed 
model we show below that the lobe size of an FRI-type source is 
not well defined, while its observed radio luminosity is likely to be 
constant for most of its lifetime. 



In our model for sources with laminar jets inflating FRII-type 
lobes, we assume that the jet is ballistic when it emerges from the 
AGN and comes into pressure balance with its lobe by passing 
through a reconfinement shock (see Canto & Raga 1991 , [Fallel 
|1991| KA). After the external medium has replaced the lobe ma- 
terial in sources with weak jets, the initially ballistic jet still needs 
to adjust to the pressure now exerted on it by the external medium. 
As the bulk velocity of the jet is supersonic with respect to the jet 
internal sound speed, this adjustment still causes a shock in the jet 
flow. In this reconfinement shock a fraction of the kinetic energy of 
the jet is dissipated to thermal energy. If the fraction is large, then 
the jet is even more susceptible to turbulent disruption downstream 
of the shock. For a fixed pressure in the external gas, weaker jets 
are more prone to disruption than more powerful jets as they have 
to dissipate a larger fraction of their power to reach pressure equi- 
librium. We identify the flare point observed in FRI-type jets with 
the reconfinement shock (Laing & Bridle 2002a}. 

More powerful jets may avoid turbulent disruption, even when 
in contact with the dense, external medium. However, usually the 
lobes inflated by these more powerful jets do not come into pres- 
sure equilibrium with the external gas in the core while the jet flows 
are still young. The jets can then inflate a large lobe, but if (3 < 2, 
then the lobe pressure will eventually fall below the pressure in the 
ambient medium and pinching of the lobes at their centres can oc- 
cur iHoo da et al.|[l994"l >. In this case the eventual replacement of 
the lobe may lead to the formation of a second, small lobe of type 
FRII closer to the host galaxy than the previously inflated lobe. This 
may explain the Double-Double Radio Galaxies (DDRGs) first de- 
scribed by Schoenmakers et al. (20001, which always consist of 
very large (Mpc-scale) outer lobes and much smaller inner lobes. 

After passing through the reconfinement shock the bulk ve- 
locity of weak jets will be slow. The jets inflate turbulent lobes in 
which the jet material and the external medium will start mixing. 
This process gives rise to the similarity of the lobes and smoke 
from a chimney. For very slow jets the analogy may extend fur- 
ther as buoyancy forces and bulk motion in the external medium 
may play an important role in shaping the lobes. Good examples 
for the latter effect are the narrow-angle tail sources where the host 
galaxies moves through the external medium creating wake-shaped 
radio lobes (e.g. O'Dea & Owen, 1986). The turbulent distortions 
of the lobe combined with the progressive mixing of the jet mate- 
rial with the external gas, make it difficult to define the lobe size 
using radio observations. Furthermore, the lobes of FRI-type ob- 
jects fade from the flare points outwards and only very deep radio 
observations can map a significant fraction of their total extent (e.g. 
|Blundell & Rawlings|[200T) . 

The flare point is usually the brightest feature of an FRI-type 
source. The relativistic electrons emitting synchrotron radiation are 
mostly accelerated here. If we assume that a fixed fraction of the 
pressure inside the jet is provided by the magnetic field and the 
synchrotron emitting relativistic electrons, then the radio luminos- 
ity of the lobe with volume V is roughly given by 

LvtI^V. (16) 

As discussed above, the size of the lobe and therefore its volume 
cannot be easily established from observations. Concentrating on 
the adiabatic inflation of the lobe immediately downstream of the 
flare point where the external pressure is roughly constant, we get 
from energy conservation 

dU=(rj-l)podV = - Pl) dV + Qdt, (17) 
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where Tj is the adiabatic index of the jet material. After a time t, 
the lobe has a volume 



Qt 



(18) 



3/4 

and the radio luminosity of the lobe is L v oc p ' Qt. This result 
neglects any energy losses of the electrons and does not take into 
account the drop in pressure in the environment once the lobe ex- 
tends beyond the core radius a. Both effects are responsible for the 
diminishing brightness of the lobe at increasing distances from the 
flare point. In practice, for a given surface brightness sensitivity 
these outer structures will contribute a decreasing fraction of the 
overall luminosity for increasing distance from the flare point. 

In our simplified model the dimming of the outer regions of 
the lobe puts an upper limit fi; m on f, the time during which the 
observable lobe was inflated. If we observe the lobe at time /, then 
those parts inflated at times < t — t\\ m will not contribute to the lu- 
minosity. The limit t\[ m depends on the pressure and the strength of 
the magnetic field in the lobe. Hence it depends on the pressure dis- 
tribution of the external medium, but not on the jet power. Once the 
lobe has reached a volume V = Qt\{ m / {TjPo) > its observable radio 
luminosity should remain roughly constant in this simple picture. 
Note that its total luminosity may in principle continue to slowly 
increase, but this is not detectable in an observation with a given 
surface brightness limit. Note also that the observable luminosity is 
linearly proportional to the jet power Q. 

Clearly the ideas outlined above do not consitute a proper 
model of sources with an FRI-type morphology. For example, we 
neglect the possibility of particle acceleration in the turbulent lobes. 
Nevertheless, based on the discussion above, in the following we 
will use the reasonable results that the observable radio luminosity 
is roughly constant and linearly proportional to the jet power. 



4 THE RADIO LUMINOSITY FUNCTION (RLF) 

The RLF of extragalactic radio sources is commonly parameterized 
as a broken power-law of the form 



n(Ly) =;j 



(hid) 



(19) 



where ij ow and i^igh are the power-law exponents at low and high 
luminosities respectively. The transition between these two regimes 
occurs at the break luminosity L„ . For convenience, the RLF is usu- 
ally expressed as number density of sources per unit logarithm of 
the radio luminosity. We follow this convention. 

The break luminosity is close to the luminosity separating the 
two FR classes with FRI-type sources dominating below the break 
and FRII types mainly found at luminosities above L* jFanaroff] 
|& Riley| |1974| >. Note however that the luminosity separating the 
two classes appears to depend on the optical luminosity of the host 
galaxy ( |Ledlow & Owe n 1996). For now we assume that the break 
in the RLF at L* is fixed and that the slopes of the RLF below 
and above the break, ij ow and ^high> are determined only by FRI- 
type and FRII-type sources respectively. We will return to the more 
complicated, observed separation of the two classes in terms of ra- 
dio and optical luminosity further on. 



4.1 The low luminosity end of the RLF 

We first concentrate on the low luminosity end of the RLF where 
n oc Ly . Above we suggest that sources with an FRI-type mor- 
phology evolve out of small FRII structures for sources with weak 
jets. Let S(Ly) be the rate per unit volume and per unit luminosity at 
which sources with FRI-type morphology and luminosity L v appear 
in the universe. The time it takes sources described by our fiducial 
model parameters, but with a jet power equal to <2min to evolve 
through their initial FRII stage is 5 x 10 6 years. Most sources that 
disrupt will have jets less powerful than this and therefore dis- 
rupt earlier. This is short compared to their total lifetime of about 
10 8 years. We therefore expect n(L v )L^ 1 oc S(L V ). The factor Z^T 1 
arises from the definition of n(Ly) as the density of sources per unit 
logarithm of Ly. 

We argued in section [372] that the radio luminosity of sources 
after they develop an FRI-type morphology is constant and linearly 
proportional to the jet power Q. We can therefore relate S(Ly) to 
the rate per unit volume and per unit jet power at which AGN pro- 
ducing jets with a given Q are appearing in the universe, S'(Q), as 
S'{Q) 01 S(Ly). Therefore the 'brith rate' of jets from AGN with a 
jet power Q must be S'(Q) °= g -1 ^ 1 ™- 

There is no a priori reason to assume that this form of S'(Q) 
also applies at the higher jet powers of the sources dominating the 
RLF above However, in the following we will show that this 
assumption leads to a self-consistent explanation of the shape of 
the RLF. For further discussion of this point see sections [6. 2. 1| and 
[6321 



4.2 The high luminosity end of the RLF 

At luminosities Ly > L„ the RLF is dominated by powerful sources 
with FRII-type morphologies. For these we have to take into ac- 
count the luminosity evolution of their lobes. Consider for now 
sources with a luminosity evolution described by L v »= /~ 5 , where 
8 will take different values according to which regime the source 
is currently in. The evolution of the RLF above is then given by 
(e.g. |Cavaliere et al.[[T97l} 

dN d ( , 9L V N 



dt dLy 



N- 



dt 



= S(Ly), 



(20) 



where N is the number density of radio sources per unit luminosity 
and so n = ln(W)L v N. 

We assume that the RLF changes only on cosmological 
timescales, i.e. on timescales much longer than the lifetime of in- 
dividual sources. This allows us to set dN/dt = 0. The solution of 
the differential equation is then given by 



Ly ^ ^ J S (^peak) dLp ea k. 



(21) 



We discuss the various terms in the following. 

Sources retaining an FRII-type morphology initially expand 
in the core region of the ambient gas distribution, D <a, where the 
gas density is roughly constant. The time taken to traverse this core 
region will be short compared to the total lifetime of the source. 
Therefore we simplify the construction of the RLF by assuming 
that sources only start to contribute to the high luminosity end of the 
RLF once their lobes extend beyond the core region, i.e. when D = 
a. The discussion in section [23] implies that the radio luminosity 
in this phase of the source evolution tends to a constant L pea k and 
Lpesk x Q- I n some sources synchrotron losses may not become 
significant during this phase. For these sources equation {3} implies 
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that L pea k oc g 7 / 6 for D = a, very close to the result for a constant 
luminosity. In the following we use L pea k oc Q for all sources. 

Once the lobes extend beyond the core region, their lu- 
minosity evolution is described by equation {3} and so 8 = 
(7p- 8) / [4(5 — (3)]. As long as [3 > 8/7, the luminosity of the 
source always decreases. Hence the luminosity L pea k at D = a is 
indeed the maximum luminosity of the source at any time during 
its life. The integration of the source function in equation \2\ \ sums 
up the contribution to the RLF of all sources with /jpeak 

Most sources retaining an FRII-type morphology will reach 
within their lifetime the regime where inverse Compton losses of 
the electrons in the lobe become important. In this regime the lumi- 
nosity evolution accelerates as 8 changes to 8' = (4 + P) / (5 — P) 
(see equation^}. We will show below that sources in the inverse 
Compton regime do not make a significant contribution to the RLF. 
Hence the upper limit of the integration is set by L max = L pea k, 
where L pea k is taken for those sources which enter the inverse 
Compton regime when their luminosity is equal to Ly. 

Using the results from above that L pea k oc Q and 5 (L pea k) ° c 
S 1 (2) x Q -1 ' bw we can solve the integral, 



-Clo 



.-(s+ivs-Uw 



(22) 



where for the last step we assume L max 2> Ly. Lpeak is the maximum 
luminosity of a source at any point during its lifetime at a time when 
D ~ a and L max is equal to Lp ea k for the sources with the most pow- 
erful jets contributing to the RLF at Ly. Their luminosity decreases 
from Lmax proportional to Z) - "( 5- P)/ 3 . These powerful sources will 
enter the regime dominated by inverse Compton losses only once 
their lobes have grown in size well beyond a. At this point their 
luminosity has decreased significantly compared to L max . 

In principle there is a second contribution to N from sources 
already in the regime dominated by inverse Compton losses. It is 
easy to see that this additional contribution has the form 



_ I ~ w 
max ^i; m 



(23) 



where L\i m > L max is now the maximum luminosity of any source 
at any time. Comparing with the term describing the sources in the 
adiabatic regime, we can neglect this contribution to the RLF at L v 
as long 3.S i^max >L V . 

With the definitions of n and N we now have 



,-^higi, , , l-(8+l)/8-Q 0W 

(lOtLy OC L V N OC Ly , 



4.3 The break in the RLF, U 

Our model can also provide an estimate for the location of the break 
in the RLF. Consider the FRII-type sources with the weakest jets. 
Their lobes just escape the core region without suffering turbulent 
disruption. Their radio luminosity then follows the relation given 
in equation {3} until the energy density inside their lobes becomes 
comparable to the energy density of the CMB, mcmb- The source 
then enters the inverse Compton dominated regime and its lumi- 
nosity evolution steepens. From equation jA9l it follows that at this 
time p = 7«cmb /12 and using equation l |A4| the source age is 

/ 12/, \^W + V (p _ 5)/3 



?IC V 7mcmb J 



s3/(4+P) ( 2 _p)/(4+p) 



X [pa?) 

The luminosity of the lobe at tic is found from equation {3} as 



(26) 



Lie 



3fL .(P- 5 )/ 3 /7 MC mbY 7F8)/[4(4+P)1 



3 + oi 1 



\ 12 / 



« ip« p ) 3 " 4+ %r +B 



(27) 



Note that for P ~ 2, Lie is roughly linearly proportional to j2min 
while tic is independent of the jet power. All sources in comparable 
environments enter the inverse Compton dominated regime at about 
the same age, independent of their jet power. However, sources with 
Q = 2min hav e the lowest radio luminosity and remain the least 
luminous objects at t > fie- 

After entering the inverse Compton dominated regime the ra- 
dio luminosity of the lobe decreases further according to equation 
(7} with Ly oc f(-*-P)/(5-P). 

Therefore, if the jet shuts down at a 
time /fl na i ~ 10 8 years, we have 



Lfinal — Lie 



( ^finaA 

V he ) 



(-4-p)/(5-p) 



(28) 



We expect the break in the RLF to occur between L nna i and Lie, i.e. 
Lfinal < L* < Lie- Also, variations in the properties of the external 
density profiles will tend to further broaden the range in luminosity 
over which the RLF changes its slope. For our fiducial model at 
v = 15 1 MHz, with Q min = 3.3 x 10 37 W derived above and f fina i = 
10 s years, we find tic = 2.8 x 10 7 years, L final = 2.3 x 10 25 WHz" 1 
and Lie = 2.9 x 10 26 WHz~' in agreement with the break of the 
RLF derived from observations (Willot tet al.||200T) . 



(24) 



or, equivalently 

_ 1 _ 4(5-p) 

Singh Slow — g — 7(3 _ 8 ' 



(25) 



Willot t et al.| ( |2001| > construct the RLF from observed complete 
samples of extragalactic radio sources. They find 2j ow ~ 0.6 and 
Chigh ~ 2.4 which is consistent with P ~ 2. The mean value of P 
for the exponent of the external density distributions of the gaseous 
haloes of elliptical galaxies, galaxy groups and clusters in the ex- 
tensive sample of Fukazawa et al. (2004) is 1.6 which would im- 
ply a larger difference i^high — Clow However, they also find that P 
is correlated with the temperature of the X-ray emitting gas in the 
sense that hotter systems have steeper density profiles. The RLF we 
study here is composed of the most luminous radio sources in the 
universe, which therefore are preferentially located in hotter and 
denser systems with steeper density profiles. 



5 THE SIZE DISTRIBUTION OF FRII-TYPE SOURCES 

The size of lobes with an FRII-type morphology increases with life- 
time of the jet flow. However, factors like the jet power and the den- 
sity distribution of the external medium influence the rate of growth 
(see equation |A2| l. Studying the size distribution of FRII-type lobes 
in flux-limited samples in detail is further complicated by the possi- 
ble cosmological evolution of these quantities as well as the overall 
lifetime of the jets and the selection effects introduced by the flux 
limit in the cosmological context (e.g. Neeser et al.| 1995 , Kaiser 



& Alexander 1999 , Blundell & Rawlings 1999 , Wan et al. 



2000; 



Barai & Wiita| |2006| 12007^ 7 Here we restrict the discussion to a 



very general study of the size distribution predicted by our model 
without taking into account cosmological effects. 

|Q'Dea & Baum| ( |l997) point out that compact sources with 
sizes below a few kpc are so numerous in observed samples, that 
not all of them can evolve into sources with lobes tens and hundreds 
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of kpc across. The number of sources per unit logarithm of lobe 
length, n(D), is constant below a few kpc, while n(D) oc D 0A for 
larger D. These findings can be explained with an intermittently 
active jet flow (Reynolds & Begelman 1997 1, but here we explain 
the same result with our model including a continuous jet flow. Our 
approach is similar to that of Alexander (2000 ) who also invokes 
jet disruption of small sources to explain the large number of small 
compared to large sources. 



5.1 Large sources 

We can use our model to predict N(D), the number of sources per 
unit lobe length, using a formalism analogous to the RLE As be- 
fore, n(D) is related to N(D) as n(D) = DN(D). Neglecting the 
possible cosmological evolution of N(D), we can write (compare 
to equation|20]for the RLF) 



(29) 



where H(D) is a sink function. As sources grow, their radio lumi- 
nosity usually decreases and so at some point they drop below the 
flux limit of the sample used to construct N(D) from observations. 
This effect is represented by H(D). 

The length of the lobe as a function of time is given by equa- 
tion i A2l as D oc f 3 /' 5 ~P'. The solution for the size distribution is 
then 



N(D) oc Z)( 2 -P)/ 3 [ Dmm ME/) dD 1 . 

■ID 



(30) 



The integration of the sink function simply sums up all sources that 
have dropped below the flux limit of the sample when their lobes 
had lengths shorter than D. Neglecting cosmological effects, the 
flux limit of the sample translates into a limit on the source lumi- 
nosity. Sources with comparatively weak jets will drop below this 
limiting luminosity while they are in the adiabatic regime. Equa- 
tion ^ then implies that the size of the lobe at the point the source 
drops out of the sample is D °= g 14 /( 7 P~ 8 ). To solve the integral, we 
note that H(D) dD °*S'{Q) dQ and so 



N(D) oc d(2-P)/3-&w(7P-8)/14 i 



(31) 



With P ~ 2 and ^i ow ~ 0.6, our model then predicts n(D) oc D 01 , 
somewhat steeper than the observations suggest. For P = 1.6, the 
mean of the |Fukazawa et al.| ( |2004| > sample, we find n(D) oc 
D l , considerably steeper than the observational value. However, 
sources with more powerful jets will enter the inverse Compton 
dominated regime before dropping out of the sample. In this case 
from equation ([7} we have D oc g( 4 +P)/ 5 and so 



Ar(D)ocD( 2 -P)/ 3 -^w( 4 +P)/ 5 , 



(32) 



or n(D) oc D - 3 for (3 = 2 and n(D) oc D 05 for P = 1 .6. For a sample 
composed of sources from both regimes, the model therefore results 
in a slope close to that derived from observations. 

For very large lobe sizes the fact that jets will stop to supply 
energy to the lobes at / = ffj na ] will become relevant. If fg na i is simi- 
lar for all sources, then equation l |A2[ l relates the final length of the 
lobe to the jet power as D oc g'/ts^W. We then find n(D) oc D~ 08 . 
The data point containing the largest sources in |Q'Dea & Baurn] 
{1997} may indicate this predicted drop in the size distribution. 



5.2 Small sources 

At small lobe sizes we expect to find many objects that currently 
have an FRII-type morphology, which will develop turbulent jet 
flows at a later stage. From the discussion of the transition from 
FRII to FRI morphology in section |3.1| it is clear that for most 
sources the transition will occur while they are contained within 
the core radius of the external density distribution, a. The density 
distribution quickly steepens outside a and as discussed above, it 
becomes easier for sources to avoid disruption in more stratified 
environments. From these considerations it follows that we expect 
the change in the slope of n(D) to occur around D — a. Given that 
n(D) as determined from observations has its break at a few kpc, 
this is consistent with the results of X-ray observations of the gas 
in elliptical galaxies (e.g. Fukazawa et al. 2 004|>. 

The density distribution given by equation \\2\ can be approx- 
imated by power-laws with an exponent changing as a function of 
distance from the centre. Hence we expect that the jets in most 
sources become turbulent when their environment is described by 
a power-law with an exponent p tul -t, < p. Furthermore, at r = a the 
exponent of the approximating power-law is P/2. If sources mostly 
disrupt inside a, we can further constrain the relevant power-law ex- 
ponent to P tur [, < P/2. Also, during phase 1 of the source evolution 
the radio luminosity increases until reaching a constant maximum. 
The sample of |Q'Dea & Baum] fl997 ) contains only very luminous 
objects. Hence many objects with weak jets disrupt before reach- 
ing a luminosity above the sample limit. Only sources with power- 
ful jets are included in the sample and this happens only once they 
have evolved long enough so that their luminosity exceeds the sam- 
ple limit. If the jets in these objects subsequently disrupt, they must 
do so at a time when the lobes have grown to sizes comparable to 
a. Hence we expect p tur b ~ P/2. With this constraint in mind we 
can construct n(D) from our model for small sources. 

We can again use equation 1 29 1, but this time H(D) describes 
sources whose jets are disrupted and which develop an FRI-type 
morphology. The condition for disruption is that the pressure in- 
side the lobe becomes comparable to that in the source environ- 
ment. From equation ( |A4[ > we then find that the lobe size at which 
disruption occurs is given by D oc g 2 /( 4 +Pi"*>. The solution of the 
differential equation is then 



N(D) oc £)(2-P.urb)/3-Uw(4+|3„ 11 t)/2^ 



(33) 



For n(D) ~ constant in this regime, we require p tul -b ~ 0.7, which 
is indeed smaller than P/2 ~ 1. For P = 2, the density distribution 
is approximated by a power-law with exponent P tul -t, = 0.7 at a dis- 
tance of r ~ 0.7a from the centre. 

In the framework of the model the flat shape of n(D) for small 
lobes is caused by many small objects developing an FRI-type 
morphology before growing to large sizes. In section [43] we have 
shown that the break in the RLF is located around 10 2t> WHz~' 
at 151 MHz. Therefore not many sources with an FRI-type mor- 
phology can contribute to the RLF at luminosities in excess of this 
threshold. However, the small sources discussed in O'Dea & Baum| 
|l997| have luminosities in the range 1 26 W Hz 1 to 10 2S WHz -1 
at 5 GHz which, for a spectral index of a = 0.75, roughly translates 
to 10 27 WHz -1 to 10 29 WHz" 1 at 151 MHz. While some sources 
with an FRI-type morphology and a radi o luminosity in excess of 
5 x 10 26 WHz 



at 151 MHz exist (e.g. Ledlow & Owen 



1996 



|Hardcastle et al.| |2007[ l, they are rare. Therefore our model re- 
quires that the radio luminosity of most sources developing FRI- 
type lobes must drop by at least two orders of magnitude once 
their jets become turbulent. Without a detailed model for the FR 
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transition it is impossible to show whether this is realistic or not. 
However, we note that this requirement may make our explanation 
for the small lobe end of n(D) problematic. Alexander (2000) also 



suggests that the luminosity of sources with disrupted jets decreases 
dramatically. However, in that model sources with disrupted jets are 
not expected to evolve into FRI-type objects, but fade so strongly 
that they become virtually undetectable. 



6 CONNECTION TO OTHER WAVELENGTHS AND 
MICROQUASARS 

6.1 Host galaxies 

In the previous sections we concentrated exclusively on the radio 
properties of radio galaxies and radio-loud quasars. However, it is 
well known that the two FR classes also differ in respects other 
than their radio properties. Ledlow & Owen ( 1996) show that the 
dividing line in terms of radio luminosity between the FR classes 
is a function of the B-band luminosity of the host galaxy. In our 
model the dividing luminosity between the FR classes is identified 
with the break in the RLF, L*. From observations we therefore ex- 
pect that L* °= Lg 8 , but note that [ 



Best et al. 



(2005) do not find a 



dependence of L* on Lb . 

In our model L* is determined by the luminosity of sources 
with the weakest jets entering the inverse Compton regime, Lie, an d 
their luminosity when their jets cease to supply energy to the lobes, 

we find Lie °=p 3/2 a 3 L 3 / 2 



15 i and i 27 



ifinal- Combining equations i 
for P = 2. Lfi na i shares the same parameter dependence and hence 
. The same parameters describing the density dis- 
tribution in the source environment also determine the X-ray lu- 
minosity of the ambient gas. Bremsstrahlung emission depends 
strongly on the density of the radiating gas and so we expect the 
total X-ray luminosity of the source environment to be dominated 
by the contribution of the core region. In this case Lx <* p 2 a 3 \/L. 
At least in elliptical galaxies the X-ray luminosity is correlated with 



1/2 i 

the optical luminosity of the galaxy as Lb oc LJ (O' Sullivan et al 



2001 1 and so we expect Lb « pa 3 / 2 ?" 1 / 4 . The temperature of the X 



ray emitting gas does not vary strongly between individual objects 
(e.g.|Fukazawa et al. , 2004), but the core radius and the central den- 
sity can be different. If the central density p was fixed for all galax- 
ies hosting powerful jets and only the core radius is varying, then 
we would expect L* °c Lg . In the case of a fixed core radius a and a 

3 12 

variation in the central density, the expectation would be L„ L fi . 
The observed correlation lies between these possibilities and our 
model is thus consistent with observations (see Gopal-Kris hna &| 
|Wiita|[2001| for a similar argument). 



6.2 Nuclear emission 

Differences between the FR classes are also found in the optical 
emission of the jet-producing AGNs themselves. Most FRII-type 
sources show strong optical line emission while FRI-types only 
possess weaker lines (e.g. |Hine & Longair| |1979) . FRI-type ob- 
jects and some FRII-types with weak emission lines are collectively 
referred to as Low Excitation Radio Galaxies (LERG) while the 
objects with strong emission lines are called High Excitation Ra- 
dio Galaxies (HERG). We emphasise that the FR classes cannot be 
identified directly with the LERG and HERG groups. There are a 
number of objects with weak emission lines, but FRII-type radio 
morphologies (Barthel 1989}, and a few objects showing strong 



emission lines combined with an FRI-type morphology (e.g. |BlurT| 
|dell & Rawlmgsl[200l) . 

The optical continuum emission of LERGs is correlated with 
the luminosity of their radio cores and is consistent with origi- 
nating in the inner jet flow rather than the AGN itself (Chiaberge 
|et al.| [2000). In HERGs the optical continuum luminosity exceeds 
the value predicted by this correlation for a given radio luminosity 
implying a contribution by the AGN itself. These results can be in- 
terpreted as evidence for a radiatively inefficient accretion flow in 
LERGs compared to more efficient flows in HERGs ( Falc ke~et al.| 
|2004[ ). This interpretation is also consistent with the observed dif- 
ferences between the two classes in X-rays (HardcastIFetaL]|2006]( 
and infrared emission (Miiller et al. , 2 004| >. 

In our model, we can combine the emission properties of the 
AGN with the radio morphology classification, if we assume that 
the radiatively inefficient accretion mode is associated with the pro- 
duction of weak jets. The main ingredient determining whether or 
not a jet disrupts and forms a lobe with an FRI-type morphology is 
the jet power. Hence with the above assumption we would expect 
AGN with radiatively inefficient accretion, i.e LERGs, to produce 
only weak jets which are more easily disrupted and give rise to ra- 
dio lobes with an FRI-type morphology. At the same time, HERGs 
should produce powerful jets which retain their FRII-type morphol- 
ogy. However, as we have seen above, the transition in radio mor- 
phology of AGNs is also influenced by the properties of the gaseous 
halo of the host galaxy. Unless the jet power associated with the 
change from radiatively efficient to inefficient accretion is for some 
unknown reason fine-tuned to exactly coincide with the jet power 
connected to the change in radio morphology, we would expect to 
observe a number of hybrid objects. We have already mentioned 
that some LERGs possess FRII-type radio lobes and it is interesting 
to note that their radio luminosities are close to L» l |Barthel|[T9"89"l >. 
Furthermore, their precarious position between laminar jet flow and 
turbulent disruption of the jets may give rise to their prominent jets 
and weak radio hotspots (Hardcastle et al. , 1998 1. Finally, HERGs 
with FRI-type radio lobes appear to be rare (e.g.[Blundell & Rawl-| 
|ings||200T| >, suggesting that the jet power associated with the switch 
from radiative inefficient to efficient accretion is lower than the 
Gmin f° r typical radio source environments. 



6.2.1 Fuelling the AGN 

Recently Hardcastle et al. (2007 1 suggested that the accretion flows 
in LERGs are fundamentally different from those in HERGs. In this 
scenario LERGs are fuelled by comparatively hot, but radiatively 
inefficient gas which is cooling out of the galactic atmosphere or 
the surrounding group or cluster environment. The fuelling by hot 
gas cooling out of the galaxy atmosphere allows for a feedback loop 
regulating the cooling of gas and therefore limiting the growth of 
the host galaxy ( Best et al. 2006 1. In contrast, HERGs are powered 
by the accretion of cold gas, most likely acquired in a merger of the 
host galaxy with a gas-rich companion. The latter accretion mode 
releases more energy in the form of radiation and, in our model, 
produces more powerful jets than the former. 

The idea of different accretion modes giving rise to differ- 
ent radio morphologies is interesting, but it also creates a problem 
for our construction of the RLF. The mechanisms for fuelling of 
sources are quite different in the two accretion modes. The switch 
between the modes occurs near to L„ and so it is not clear whether 
the birth rate of jets, S'(Q), determined from the low luminosity 
end of the RLF in section [4~T| can be extrapolated to the high lu- 
minosity end as we have done in section [4~2] The overlap in radio 
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luminosity of LERGs and HERGs around L„ makes a clear separa- 
tion of S'{Q) for the two classes very difficult. However, the model 
is self-consistent and produces the break in the RLF without invok- 
ing a difference in the birth functions of LERGs and HERGs. 

6.2.2 Relation to microquasars 

An alternative to the above scenario is provided by analogy to mi- 
croquasars. These objects show radiatively efficient and inefficient 
accretion states without the need for a difference in the fuel source. 
The transitions are moderated by the accretion rate through the in- 
ner parts of the disc. Instabilities within the disc can lead to large 
changes of this parameter (e.g. Frank et al. , 1992). For high accre- 
tion rates (high-soft state) the discs are luminous, but formation of 
jets appears to be suppressed, while for small accretion rates the 
discs are radiatively inefficient (low-hard state) and comparatively 
weak jets are produced ( Fe nder et al.| |2004). During the transition 
from the low-hard to the high-soft state, powerful, but short lived 
jet ejections occur. The LERGs may be identified with the low-hard 
state of microquasars as they also show radiatively inefficient accre- 
tion discs and weak jets (Kording et al. 2006 ). In this scheme, the 
HERGs must be identified with the transition between the states as 
the formation of powerful jets and a luminous accretion disc is re- 
quired for them. Radio-quiet quasars may then be the counterparts 
to objects in the high-soft state. 

Is the state transition slow enough in AGN to explain the life- 
time of jet flows? It is difficult to measure the duration of the state 
transition in microquasars, but from simultaneous X-ray and radio 
observations the relevant timescales appear to be of order several 
hours to a few days (e.g. |Fender et al.| |1999| >. This must be the 
timescale for the disc instability causing the state transition. It is 
by no means clear which physical mechanism is responsible for the 
disc instability. However, it is likely that the relevant process is con- 
nected with the viscous timescale within the accretion disc which 
roughly scales as the mass of the central black hole to the power 
3/2 ( |Frank et al.| |1992[ >. The mass of the black hole in a typical 
microquasar is lOM©, while that of a black hole at the centre of 
an AGN is about 10 8 Mq. Hence we would expect the duration of 
the state transition in an AGN to be of order 10 8 years, compara- 
ble to the maximum lifetime of jet flows in powerful radio galaxies 
l |Alexander & Leahy|[T987) . 

If this scenario is correct, then the difference in the radiative 
properties between LERGs and HERGs, and ultimately that be- 
tween the FR classes, is not caused by differences in the fuelling 
mechanism. In this case, weak and powerful jets can originate in 
the same system at different times. The extrapolation of the birth 
rate of jets S'(Q) to higher jet powers is then less problematic. Also, 
black holes of a given mass can give rise to jets with a range of jet 
powers, translating to different radio luminosities of their lobes, at 
different points during their accretion history. This predicted statis- 
tical independence of the radio luminosity and the black hole mass 
is consistent with observations (Best et al. 2005 1. 



7 SUMMARY 

In this paper we use an analytical model for the dynamics (KA) and 
the luminosity evolution (KTJA) of individual radio galaxies and 
radio-loud quasars with FRII-type morphologies to explain some 
of the observed properties of the entire source population. We show 
that the radio luminosity of the lobes evolves through phases gov- 
erned by the dominant energy loss mechanism of the radiating, rela- 



tivistic electrons. While the lobe expands in the quasi-constant den- 
sity core close to the centre of the host galaxy, the luminosity first 
rises and then levels out to a constant value as synchrotron losses 
become dominant. At the same time the spectrum steepens, consis- 
tent with the properties of the small CSS sources. As the lobe grows 
larger than the core radius, the external density decreases and syn- 
chrotron losses become less important. Consequently, the spectrum 
flattens again. 

We complement the model for FRII-type objects by presenting 
a very simplified model for FRI-type lobes. In this approach, the ra- 
dio luminosity of the FRI-type lobes is constant for most of their 
lifetime and linearly proportional to the jet power. We also pos- 
tulate that all radio galaxies and radio-loud quasars start their life 
with an FRII-type morphology. However, the lobes of objects with 
weaker jets come into pressure equilibrium with the surrounding 
gas before the lobe extends beyond the core region. Starting close 
to the centre of the source, the lobes are then buoyantly replaced by 
the denser ambient medium which will eventually reach the jet flow 
itself. Without its protective lobe a weak jet can turbulently disrupt 
and develop an FRI-type lobe. Somewhat more powerful jets may 
be able to avoid disruption and to produce another pair of FRII-type 
lobes. The DDRGs show this kind of structure. 

With the combined model for both FR classes we can self- 
consistently explain the broken power law shape of the RLF. The 
predicted break luminosity is close to that derived from observed 
complete samples. The RLF below the break luminosity is assumed 
to be dominated by FRI-type objects. The slope of the RLF in 
this region gives the rate at which AGNs producing jets with jet 
power Q appear in the universe. This birth function is proportional 
to Q . Above the break luminosity, FRII-type objects dominate 
and the observed slope of the RLF can be successfully reproduced 
from the shape of the birth function and the known luminosity evo- 
lution. 

We also recover the observed distribution of lobe sizes for 
FRII-type objects. The surprisingly flat size distribution at small 
sizes is caused by the increasing number of sources developing 
FRI-type lobes and thereby 'dropping out' of the distribution. How- 
ever, this result depends critically on a significant decrease of the 
radio luminosity of those sources transforming into FRI-types at 
the point in time their jets disrupt. 

We explain the dependence of the radio luminosity delineating 
the FR classes on the optical luminosity of the host galaxy. More 
luminous hosts are located in denser environments in which more 
powerful jets can be disrupted. More powerful jets would normally 
lead to more luminous lobes and so the demarcation line between 
the FR classes in terms of radio luminosity increases with increas- 
ing optical luminosity of the host. 

Finally we discuss the connection between the recently pro- 
posed difference in the accretion modes of LERGs and HERGs and 
their radio morphologies in the context of our model. Our findings 
are consistent with the idea of radiatively inefficient accretion of hot 
gas in LERGs and luminous cold gas accretion in HERGs, if radia- 
tive efficiency is correlated with jet power. However, differences in 
the fuelling mechanism between the two classes may pose a prob- 
lem for our construction of the RLF as we assumed the same birth 
function for all jet-producing AGN. The problem can be avoided, 
if the radiative efficiency of the accretion process is regulated by 
instabilities in the accretion disc itself, independent of the fuelling 
mechanism. This latter idea is inspired by the observed situation 
in microquasars. It may provide yet another link between accret- 
ing and jet-producing black holes spanning orders of magnitude in 
black hole mass. 
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APPENDIX A: RE- WORKING THE KA AND KDA 
MODELS 

In the following we derive expressions for the dynamical evolution 
and radio synchrotron emission properties of FRII-type radio-loud 
AGN. Unless stated otherwise, these expressions are identical to 
the relations given in the papers Kais er & Alexander| < fT997| KA) 
and |Kaiser et~aT1 ( [T99T . KDA). Our reason for re-deriving them is 
to demonstrate how they depend in detail on the most fundamental 
properties of the jets and their gaseous environments. In particular 
we concentrate on the jet power, Q, the length of the lobe along the 
jet axis, D, the age of the jet flow, t, the pressure inside the lobe, 
p, the monochromatic luminosity density, L v , at a given observing 
frequency v and the 'density parameter', pa~P (see below for an 
explanation). 

In what follows we necessarily need to frequently reference 
material in the KA and KDA papers. We abbreviate references in 
the following way: KA4 refers to equation 4 in KA, while KDA2 
refers to equation 2 in KDA. 



Al Source dynamics 

In our discussion we concentrate on a single lobe inflated by one 
of the jets in a radio-loud AGN. The jet is assumed to come into 
pressure equilibrium with its lobe through a reconfinement shock. 
After this, the jet flow remains laminar until it enters the hot spot 
region where its kinetic energy is dissipated. We also assume that 
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the density distribution of the medium external to the lobe is given 
by a power-law distribution according to 



Px = P 



(Al) 



where p and a are constants and r measures the distance from the 
AGN assumed to be located at the centre of the distribution. The 
power-law form of the density distribution implies that p and a are 
not independent parameters and that the model can only depend on 
their combination pa$, which we refer to as the density parameter. 

With these assumptions, the length of the lobe is given as 
(KA4 and KA5) 



D = c\ 



1/C5-P) 



-3/(5-P) 



The constant c\ can be written as 

(A 4 (r, + l)(T,-l)(5-P) 3 



CI 



18ti 9[r, + (r,-l)A 2 /2] - 4-p 



V(5-P) 



(A2) 



(A3) 



The ratio of the lobe length and its radius, A, stays constant in our 
self-similar model. Here, we assumed a cylindrical geometry of the 
lobe. Note that KA uses the aspect ratio of the lobe which is 
denoted R in KDA. Our choice for A implies A = 2Rj = 2R. T/ and 
r x are the adiabatic indices of the lobe material and the external 
gas, respectively. Note that this expression for c\ differs slightly 
from KA25 due to an error in KA first pointed out by Kom issarov| 
|& Falle|p998) . 

The pressure inside the lobe is given as a function of the source 
age / by KA20. Using equation \h2\ and following the discussion 
on the expansion of the lobe leading up to KA38, we find 



P =f p (paty /3 Q yw-*-v/\ 



with 



fp 



18c 



2(5-P)/3 
1 



(r,+i)(5-p)*A2 



(A4) 



(A5) 



Note here that by substituting equation |A2| into equation {A4|, we 
can express the pressure in the lobe as a function of the source age, 



,M'f" b b(h»/<h>>,(-«-i»/<!-i». 



(A6) 



This shows that if P ~ 2, then at a given point in time the pressure 
inside the lobe is the same in all sources located in environments of 
similar density, regardless of the jet power. 

A2 Synchrotron emission 

We now turn to the processes giving rise to the radio synchrotron 
emission of the lobe. We have to take into account the energy losses 
of the relativistic electrons inside the lobe. These depend on the 
time the electrons have spent inside the lobe and so depend not 
only on the age of the jet flow, t, but also on the time t\ < t at 
which they were injected into the lobe. The approach of KDA is to 
split the lobe into small volume elements dV which are 'labelled' 
with their injection time f;. The evolution of the electron population 
inside these dV can then be traced individually in a self-consistent 
manner. KDA show that for adiabatic expansion, dV « t" 1 with (see 
KDA 14) 



4 + P 

r/(5-py 



(A7) 



The radiating relativistic electrons are accelerated at the end 
of the jet flow such that their total energy and the energy stored 
in the magnetic field at least initially follow the requirements of 
the minimum energy condition (e.g. Longair, 1994). The energy 
distribution of the relativistic electrons follows a power-law at time 
t[, when they are injected into the lobe, 



n(Y)dy=n (fi)Y m dY, 



(A8) 



with y the Lorentz factor of the electrons. At the time t[ the energy 
densities of the magnetic field, u^, and that of the electrons, u e , is 
then fixed to (KDA15) 



I( B (ti) 



U e {ti) = 



P(tl) 



(r,-i)(*+i)(M-i)' 

rp(k) 



(A9) 



(A10) 



(r,-i)(r+i)' 

where r = (m+ 1) /4 and k is the ratio of the energy stored in 
non-radiating particles and the sum of the energy in the magnetic 
field and the relativistic electrons. Since p(t) = p(t\) (r/fi)~ r ' fll and 
m b(0 = "B( ? i) (t /ti)~ 4a '^ , the energy density of the magnetic field 
in the lobe at time t is 

ai(r,-4/3) 



MO 



p(0 



(All) 



(T,-l)(*+l)(r+l) \ti 

The energy density and hence the strength of the magnetic field is 
not uniform throughout the lobe, if T/ ^ 4/3. Note here that we do 
not consider the case where the magnetic pressure in the lobe fol- 
lows a non-relativistic equation of state (case C in KDA). Such a 
situation would require a very special behaviour of the MHD tur- 
bulence inside the lobe constantly removing energy from the mag- 
netic field during the expansion. Hence we set Lb =4/3 wherever 
the adiabatic index of the magnetic field appears in KDA. 

The Lorentz factor of a relativistic electron emitting radiation 
predominantly at frequency v is in SI units 



Y = 



27tra e V 



(A 12) 



where the strength of the magnetic field is related to the energy 
density of the field by 

B= V^wb, (A13) 

with fjQ the permeability of the vacuum. From the discussion of the 
evolution of u-g, it follows that 

/ N 8l (4/3-r,)/4 

y=fyp u) ' (A14) 

where the constant fj can be found from equations i A9 1 and I A12 1 

as 



fr 



\ 



27im e V (T,-l)(k+l)(r+l) 



(A 15) 



The normalisation of the energy distribution of the electrons, 



equation i A8 1, is set by the energy density u e (ti) (see KDA8) and 
can be written as 

Oil"; 



«0 (ti) =fnP 



(A 16) 



where the constant f„ is related to the maximum and minimum 
Lorentz factors y max and y m j n at the time of injection by 



fn- 



(F,-l)(r+l)m eC 2 



Imin te 



—m 
max 



m — 1 



\-m \-m 
'min 'max 



,(A17) 
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for m^2 and 

fn = 



(r z -l)(r+l)m e C 



ln- 



Yi 



mm rmax 



(A 18) 



for m — 2. 

We now use the results above to express the monochromatic 
radio synchrotron luminosity in such a way to emphasise its depen- 
dence on the most important source parameters. For this we sub- 
stitute the expressions derived above into equation KDA16. While 
KDA16 gives the radio emission per unit solid angle, P v , we assume 
that the emission is isotropic and so the radio luminosity derived 
here is given by L v = 4irP v .The result is 



Ly = flQp 



(m+1)/4 f^' -v e [l-Yg(.V)] m ~ 2 



dv. 



where 



ft 



2A 2(i-r,)/r, / 3-« /fi 



3(*+l)(r+l) v ' 

with Ot the Thomson cross-section, and 

/ m \ ai / 4\ , 
e=(!+3 -Ti)ai + ^\Ti--\{3-m) 



(A 19) 



(A20) 



(A21) 



We have also changed the integration variable from /; in KDA16 to 
x = ti/t. The part of the lobe injected at the time /; corresponding 
to JCnjij, is the oldest part of the lobe still contributing to the overall 
emission at the current time t. The value of x m ; n is determined by 
the energy losses of the relativistic electrons and these are repre- 
sented by the function g (x,t). We will discuss this function and the 
determination of x m \ a next. However, there is no general analytic 
solution for the integral in equation l |A19[ >, 



4/3, this Lorentz factor Yv depends on the injection time f; and is 
therefore in general not uniform throughout the lobe. The energy 
of all electrons is decreasing monotonically. Therefore the oldest 
part of the lobe still contributing to the emission at frequency v 
is defined by the condition that electrons with a Lorentz factor Yv 
at the current time t had a Lorentz factor Yi = Ymax at the time they 
were injected into the lobe, t[. Thus the lower limit of the integration 
in equation |AT9}, x m [ n , is defined by 



Yv 



. , 1/4 ai(T,-4/3)/4 . 
■ JyP X min 



1 +-C / n 3 Ymax,?(*min,f) 



For very large Ymax this can be simplified to 



fyP 



-1/4, 



fli(4/3-r,)/4 
*min 

8 t-^min? 



(A26) 



(A27) 



but even in this case the relation can usually not be written as an 
explicit expression for x^ n . 



A3 Energy losses of the relativistic electrons 



Inside the lobe, the relativistic electrons are subject to energy losses 
due to the expansion of the lobe, the emission of synchrotron radi- 
ation and the inverse Compton scattering of CMB photons. Their 
Lorentz factor evolves according to 



dY 



"±1 
' 3 t' 



4<Tt 

3m e c 



"B + "CMB J 



(A22) 



where mcmb is the energy density of the CMB photons at the source 
redshift. We use the initial condition that the Lorentz factor of an 
electron is Yi at the time of injection into the lobe At time t the 
Lorentz factor of the same electron is y. The solution of this differ- 
ential equation is then given by (compare with KDA6 and KDA7) 



Y = 



W 3 Yi 



l+W 3 Yig(x,0' 
where 



. 4o T 

r M = ~ — t 

3m e c 



l 'B (Q ,j | "CMB 

«3 a/\ 



\-x ai ) 



for a\ 7^ 3/5, with £13 = 1 — 5a\ /3 and 04 = 1 — a\ /3, and 



. 4o T 

3m e c 



-«B {t)\nx+ -mcmb 



-x 4 / 5 ) 



(A23) 



(A24) 



(A25) 



for a\ = 3/5. The latter case occurs when T/ = 4/3 and p = or 
r"; = 5/3 andP = 1/2. The function g(x,t) is equal to t a '^a 2 (t,ti), 
with 02 the function defined in KDA7. 

At time t the Lorentz factor of the electrons emitting radiation 
at the observing frequency v is given by equation l |A12[ l. For F/ 7^ 



